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1. Introduction 

Anomalies are arguably among the most interesting phenomena to come out of studies of 
quantum matter. Their importance lies in their robustness across various length/energy 
scales - as one passes from one description of matter into another, anomaly matching ala 't 
Hooft ensures that the underlying anomalies of a theory survive in various disguises. 

While this statement is relatively better understood within the realm of effective theo- 
ries (as exemplified by the phenomenology of WZW term in particle physics and solid-state 
physics) , we have only a limited understanding of the role of anomalies in various finite tem- 
perature/finite chemical potential setups. Any progress in the phenomenology of anomalies 
is welcome especially given the important role of quantum anomalies and their associated 
transport phenomena in fields ranging from solid-state physics to cosmology. 
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To be more precise, we are concerned with the following situation : consider a quantum 
system with a continuous symmetry 1 which via Noether theorem corresponds to conserved 
Noether currents. Consider gauging this symmetry by introducing a set of external non- 
dynamical gauge fields. We will say the quantum system has an anomaly if in the presence 
of such non-dynamical gauge fields the covariant Noether currents are no more conserved. 
One can now consider instead turning on temperature/chemical potential for the Noether 
currents and ask what novel processes are characteristic of a quantum system with underlying 
anomalies. 

By now, such transport processes have been studied from various points of view - they 
are known to be constrained by thermodynamics/adiabaticity in arbitrary dimensions [1—4], 
the corresponding transport coefficients can be derived via a Kubo-like formula in 4d [5]. 
They are an established feature in various holographic fluid phases in CFTs dual to AdS3 
[6] and AdSs [7-11] where the CFT anomalies are in one to one correspondence with the 
Chern-Simons terms in the AdS bulk. The effect of Chern-Simons term for a U(l) gauge 
field in arbitrary AdS2n+i/CFT2 n was worked out in [3] - no analogous results are known for 
gravitational Chern-Simons terms in higher AdS spacetimes 2 . 

In 2d field theories, the relation between U(l) anomaly and transport is one the founda- 
tions of modern theories of Hall effects and there is an analogous relation between thermal 
transport and gravitational anomaly (see for example [14]). There are by now various ways 
in which such transport processes in 4d free fermion theories have been derived - some of 
them quite old [15-20] and others more recent [21-23]. Further, since 4d Weyl fermions are 
systems with Berry phases, these transport processes have close links to the general theory of 
Berry phases and transport [24, 25] and in particular transport in Weyl semi-metals [26, 27]. 
The transport processes linked to anomalies have also made their appearance in discussions 
about classification of topological insulators [28, 29]. 

At present, the most general set of results were derived via thermodynamic arguments 
employing adiabaticity[3, 4]. These results can be summarized as follows 3 : In a fluid the 



To simplify our discussion, we will assume this global symmetry is not spontaneously broken - in other 
words, we are interested in transport processes in normal fluid not superfluids. Though various ideas that 
we discuss in the article have their counterparts in superfluids, we believe the phenomenology of anomaly- 
induced transport in superfluids is sufficiently different to merit a separate discussion. There is by now a 
vast literature on such transport phenomena which are beset with their own subtleties. Since it would be 
too tedious/distracting to compare and contrast the effect of anomalies in the two situations, we will choose 
entirely focus on normal fluids in what follows. 

2 It would be interesting to construct and study rotating solutions in AdS2n+i with pure/mixed gravitational 
Chern-Simons terms and link the proposed modification of Wald entropy [12, 13] against the anomaly- induced 
entropy transport in the dual CFT. 

3 We will use the notations of [4] in the following. See section §^ and the appendix §|5] for a discussion of 
the basic setup and notations. 
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energy, charge and entropy transport are given by 

TT ee eu»u v + pP^ + q» anam u» + u^ anom + TZ S 

= UU^ + J anom + J diss (1-1) 
Jg EE SU^ + J^anom + ^S,diss 

where is the velocity of the fluid under consideration which obeys u^u a = —1 when 
contracted using the spacetime metric g au . Further, P^ v = g^ v + u^u v , pressure of the 
fluid is p and {e, n, s} are the energy,charge and the entropy densities respectively. We 
have denoted by {qanom, Janom, J§ anom \ the anomalous heat /charge/entropy currents and by 
{ T dLi J diss> J s,dis S } the dissipative currents. 

We are primarily interested in the anomalous currents in what follows. It is convenient 
to work with forms - let {q a nom, Janom, Js,anom} be the Hodge duals of the corresponding 
currents 4 . Using adiabaticity the following statements can be made in flat spacetime[4] 

1. All these currents are derivable from a single Gibbs current Qanom which describes the 
transport of Gibbs free energy G = E — TS — uQ in the fluid. We have the following 
thermodynamic relations 

dQ anom 



Janom 



Js 



axiom (l-^) 



dT 

Qanom = Qanom ,anom ~\~ H Janom 

2. Qanom is determined in terms of the fluid vorticity 2-form u , the rest-frame magnetic 
field 2-form B and a (n + l) th degree polynomial ^ nom \T, n] in temperature T and 
chemical potential fi (where d = In is the number of spacetime dimensions). The 
explicit expression for Qanom is given by 5 



(2w 



K> {r anom [T{2u), B + M 0)] - [^ nom [T(2a;), B + a(2u)\ 



uj=0 



15uj anom 



[T(2u),B + n(2u)]] } 



(1.3) 



A u 



3. The polynomial $anom[T, lA obeys two constraints - first, it has no term linear in T. 
Second, its value at zero temperature is completely determined by the U(l) anomaly in 
the system. 



4 Throughout this article, we will use overbars to denote Hodge duals. 

5 This is just a convenient rephrasing of the formulae/results presented in the Appendix §A of [4]. In 
particular, we have the following relation relating the functions appearing there to the functions appearing 
here 

^„ om [T(2w), B + /i(2w)] = f[B + /i(2w)] + V (2lo) 2 q[B + M (2w), Tu] 
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4. To these statements, we can add the following statement which does not follow from 
thermodynamic arguments in [4] but nevertheless seems to be true across various sys- 
tems - the T 2 coefficient in 2d and 4cZ seem to be related to gravitational anomalies (see 
for example [11, 23] for results in Ad. 

This summary encompasses everything that is known till now about these transport 
processes just from flat spacetime thermodynamics alone. But, this is unsatisfactory for 
various reasons - first of all, as we had mentioned various gravitational anomalies of the 
system show up even in the flat spacetime transport and this is mysterious even from the 
point of view of flat spacetime thermodynamics. Second, we do not know how these relations 
to gravitational anomalies generalize to higher dimensions. Further, a more microscopic 
understanding of these transport processes would clearly be useful for various reasons - for 
example, we would like to study these transport processes away from equilibrium/in the 
presence of disorder/ in the lattice analogues of continuum Weyl fermions. Our aim in this 
article is to begin addressing these questions in the simplest system exhibiting such transport - 
free fermion theory in even spacetime dimensions d = {2n — 1) + 1 with a collection of fermions 
with different chiralities Xd=2 n an d charges q under some specific £7(1) global symmetry of the 
free theory. 

After establishing the basic setup/notations in the section we proceed to study in 
detail the simplest case of 1 + Id chiral fermions in flat spacetime in section §||. We then 
translate the thermodynamic arguments of adiabaticity in flat spacetime into a more micro- 
scopic set of equations which one can take as an alternate starting point for the 2d analysis 
and which easily generalizes to higher dimensions. We propose the following intuitive picture 
for how the transport phenomena linked to anomaly arise : first of all, we propose that these 
transport processes can be thought of as arising from a certain chiral density of states and 
their spectral flow. In a continuum description, we capture this by a chiral spectral current 
Jq(x,E p ) whose time component is the chiral density of 1-particle states with charge q and 
the rest-frame energy E p and whose spatial components tell us about the flow of such states as 
the fluid flows. In section §||, we argue that adiabaticity can be simply seen as a conservation 
type equation for this chiral spectral current. 

This is of course a very well-known aspect of anomalies whereby the basic process driving 
the anomaly is the continuous injection of chiral zero modes by the external magnetic field 
(which we will assume to be small and slowly varying) into the system. We make this 
precise in section §|| by simply matching the microscopic discussion of adiabaticity to the 
thermodynamic discussion of adiabaticity. This matching gives a boundary condition for the 
conservation equation for the chiral spectral current by giving the rate at which chiral modes 
are injected into the fluid. 



These modes, once injected by the magnetic field are then convected along with the fluid 
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flow and the adiabaticity is just the statement that there is no more creation/annihilation of 
states at finite energies. 

Armed with this intuitive picture, we proceed in section §^ to solve the conservation 
equation in flat spacetime. This results in a simple expression for the chiral spectral current 

x _ X d ^ n ( qB + 2uE p \ n_1 u 

J «-^r[ — 2^ — ) A (^i)i (L5) 

where Xd=2n ^ s the chirality of the 1-particle state. This solution tells us how the states of 
different energies and charges flow hence solving once for all the spectral flow problem in an 
arbitrary fluid flow (in flat spacetime). We now notice a remarkable result - if we take the 
external field strengths to zero B — > 0, there is no more an injection of new zero-energy states 
into the fluid, but the chiral spectral flow is still non-zero even if the anomaly is turned off 
! In this case, the vorticity is sufficient to drive the chiral spectral current and this is the 
basic reason why rotational response encodes information about the anomalies in the system 
- both gravitational and non-gravitational. 

In the next section §|7], we add up the Gibbs-free energy contribution of each 1-particle 
state to get the anomaly- induced free-energy current Ganom- We find that Qanom is of the 
form given in eqn.( |1.3| ) which was derived in [4] by thermodynamic considerations. While 
this is not surprising, we find the the polynomial danom 1S derived by a very simple formula 
from the anomaly polynomial of the system. We get 

VZnom = Tanom [F (i, Pl M- -T 2 , p k>1 (JH) M- 0) (1.6) 

where V a nom is the anomaly polynomial of the system written in terms of the U(l) field 
strength F and the k th Pontryagin forms 6 of the spacetime curvature p k (9\). The above 
formula gives a simple replacement rule by which one can go from the anomaly polynomial 
to the polynomial 3 anom- 

Note that this a generalization of the observation made in [23] that T 2 coefficient in Ad 
free theories seem to be related to gravitational anomalies. Since the 4d relation continues 
to hold in strongly coupled holographic phases too [11], it is tempting to conjecture that the 
above replacement rule would continue to hold even beyond free theories. We emphasize that 
this is quite a non-trivial conjecture and what we have is a preliminary evidence that it might 
be true. We discuss this along with other further directions in our discussion section §||. We 
collect various useful results in our appendices. In the next section, we begin by explaining 
our basic setup and defining our notation - most of it being quite standard and elementary, the 
reader should feel free to skim through the section just noting various remarks on notation. 

2. The Basic setup 

The main system we will be concerned about throughout this article is a system of free 
relativistic fermions at finite temperature and chemical potential in a flat spacetime with the 



°See appendices §M for definitions of various quantities related to anomaly polynomials. 
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spacetime dimension d = 2n = (2n — 1) + 1 being even. Every particle state occurs along with 
its anti-particle state and we will call such a particle/anti-particle pair as a species. Hence, 
we will denote by ^2 species the summation where each particle/anti-particle pair contributes a 
single term to the sum. This should be distinguished from Y2f w hich denotes the summation 
where each particle (or anti-particle) state contributes separately to the sum. 

Among the large symmetry that the free theory enjoys, we will choose a specific U(l) 
symmetry for which we will turn on the chemical potential fi and put the system at a finite 
temperature T (we will use the letter /3 = 1/T to denote the inverse temperature) . The 
thermodynamic T and fi are of course defined on a local rest frame defined via a unit time-like 
vector u^. The thermodynamic potential appropriate to such a grand-canonical description 
is the Gibbs free-energy density Q = e — Ts — fin where e, s, n are the rest-frame energy 
density, entropy density and the charge density respectively. The first law takes the form 
dQ = — sdT — ndfj,. 

For a free theory, the Gibbs free-energy is obtained by simply multiparticling contribution 
from the 1-particle sector. We will denote the Gibbs free-energy contribution of a fermionic 
state with a charge q and rest frame energy E p by g q . Hence 



9 q = -\\* 



1 + e -P(Ep~w) 



The occupation of a given 1-particle fermionic state is given by the Fermi-Dirac distribution 
denoted by f q 

f = I 

The contribution of a 1-particle state to the entropy is given by the negative of the Boltzmann's 
^-function which for fermionic states takes the form 

W, = /,ln/, + (l-/,)Iii(l-/,) 

Using the standard terminology , we will call a state with slowly varying T, fi, u M with 
the local density matrix being close to the thermal density matrix as an Ideal gas. Since 
the constituents are Weyl fermions, we will call this an Ideal Weyl gas. In this article, we 
will be concerned about a specific subset of transport processes which are linked to various 
anomalies in the ideal Weyl gas. Before proceeding let us dispose of a specific technicality : 
as is well known, in the strictly non-interacting limit, all the dissipation length/time scales 
diverge which in turn means that dissipative transport coefficients like viscosity also diverge 
in this limit. While this is true, this shall not worry us too much since the effects that we are 
looking for are non-dissipative and are not plagued by such 'free-theory' infinities. While we 
expect addition of interactions/dissipation would not modify our results (this expectation is 
partially justified by various existent calculations in strongly coupled holographic phases), it 
would be nice to explicitly prove this statement. 

Having addressed that technicality, let us continue : we want to study such an ideal 
Weyl gas in presence of non-dynamical background electric and magnetic fields. Let F be the 
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field-strength 2-form, we define the rest frame electric 1-form via E^ = F^u". We can then 
do an electric-magnetic decomposition 

- [u„E v - E^uv) = B„ v (2.1) 

or in the language of forms 

F = B + u A E (2.2) 

where B is the rest-frame magnetic 2-form completely transverse to u^, i.e., B^ u u v = 0. We 
will also use the standard decomposition of the velocity gradients 

9 . , 

D^Uv = a^v + ojfja, - Uf_.au + - - - P^ (2.3) 

in terms of the shear strain rate a^, the vorticity cj„„, the acceleration and the expansion 
rate 6 of the fluid. Further P^ = g^yUpUy as before. This in particular means the exterior 
derivative of the velocity 1-form has the decomposition 

Du = 2uj-uAa (2.4) 

where uj is the vorticity 2-form. 

Before we enter the main argument of the paper, we will make some comments regarding 
our conventions for chirality in higher dimensions. Consider the 1-particle states which are 
given by the solutions of the Weyl equation of appropriate chirality (which is just the massless 
Dirac equation with the opposite chirality projected out). We will define our conventions for 
chirality now by essentially equating it to the helicity. To do this, let us divide the 2n — 1 
spatial directions into a direction x\ and n — 1 planes where the k th plane is the (xi 2k , x i 2k+1 ) 
plane where k = 1, . . . , n — 1. 

Consider first the positive frequency solutions of the Weyl equation. Further we will 
consider only the solutions with only p° > O^p 1 ^ - all other components of momentum 
being zero and the spin along k th plane being S{ 2k i 2k+1 = \ok with = ±1. Any other 
solution can be obtained by rotating this solution or by linear combinations thereof. For 
these solutions, we will assign the chirality via their helicity 

n-l 

X d=2n = sign(pi) Yl °k 

k=l 

For a Weyl equation of a particular chirality X d =2 n onr y the obeying the above equation 
are allowed. When n = 1 (when d = 1 + Id) this gives a single state with the sign of p\ fixed. 
For n > 2, we get 2 n ~ 2 states with p\ being arbitrary. 

We now turn to the negative frequency solutions which are complex conjugates of positive 
frequency anti-particle solutions. The complex conjugate of a Weyl spinor of chirality X d=2n 
is another Weyl spinor of chirality (— l) n ~ 1 Xd=2n- Hence, for corresponding to every particle 
state above we get an anti-particle state with chirality (— l) n ~ 1 x d=2 „- We will call the Weyl 
fermions with positive chirality as left fermions and those with negative chirality as right 
fermions. This summarizes the basic definitions needed for the rest of the paper. The reader 
can find a table of notation in the appendix §[D] for ready reference. 
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3. Anomaly and transport in 2d Weyl gas 

Before trying to tackle the case of a Weyl gas in higher dimensions, it is instructive to work 
out the simplest case of 1 + 1 dimensions. This is a very well-studied system and in some sense 
we will not have anything new 7 to add except for a way of looking at the standard results 
which will prepare us for the more subtle effects in higher dimensions. With this objective in 
mind, we will focus on this simple case in some detail. 

We will begin by considering a single species of a free left Weyl fermion in 1 + 1 dimensions 
with charge q. We will assign this fermion a chirality x d =i = +1- The anti-particle of this 
fermion is again a left Weyl fermion with charge —q. This follows from the general rule that 
in (2n — 1) + Id the chirality of the anti-particle is (— l) n+1 times the chirality of the particle. 

We are interested in a gas of such Weyl fermions at a finite temperature T and chemical 
potential \i. These quantities are of course defined in a center of mass frame of the ideal Weyl 
gas - let this frame be defined by a 1 + 1-d unit time-like vector which we will take it to be 
constant. It is this which in the hydrodynamic description will describe the fluid velocity. 

The question we want to address is this - what is the hydrodynamic description of such 
an ideal gas ? We will first give a naive answer to this question which will later correct. 

The conventional intuition is that this system behaves like an ideal fluid with the following 
naive constitutive relations for energy /charge/entropy currents 

^aive = ™" (3-1) 
^S.naive = SU ^ 

where the energy density e, pressure p , charge density n and entropy density s can be 
calculated from the usual statistical mechanics of an ideal fermion gas. We will calculate this 
in a moment, but before that we will argue that the above form is definitely incomplete! 

The reason is simple - a theory of a free Weyl fermion is a holomorphic 2d CFT and hence 
only the holomorphic components of the currents can be non-zero. The above relations are in 
clear contradiction with holomorphy - for one, the charge/entropy currents are time-like rather 
than null as would be predicted by holomorphy. So we are led to the surprising statement 
that conventional semi-classical intuition about the ideal gas is in direct contradiction with 
holomorphy in this simple system. Having concluded thus, let us actually calculate carefully 
what the actual constitutive relations should be. 

A left Weyl fermion field in 2d is just a single component complex field ip which obeys 
the Weyl equation (or the massless Dirac equation) 

[dt + d x ]4, = o 

This follows from the particular choice for the Gamma matrices {r*,r x } = {—ia y ,a x }. We 
can repeat the same exercise for the right Weyl fermion with x d=2 = ~~ 1 where we just flip 

7 See for example [30] for a different take on anomaly /transport in 2d fluids. 
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the sign of the d x term. In the following we will write our formulae in such a way that the 
expressions for the right Weyl fermion can be obtained by putting x d=2 = ~ 1- Hence we write 
the Weyl equation as 

[d t + Xd=2 d x ]i> = 

This equation is easily solved - solutions are just plane- waves that travel from left to right in 
the space 



- f°° —L I \ a Jp-x + fct e -v- 



P»=E p {l, Xd=2 } 



where a P and bp are the creation operators for the particle and the anti-particle respectively 
and E p is the energy in some arbitrary frame. 

Let us work in the rest frame define by from now on - so we take = {1,0}. Let e^ v 
be the completely antisymmetric tensor in 2d with e tx = 1 which implies e^ v u v = {0, 1}. The 
Weyl equation is 

and the plane-wave solutions above are 

JO 2 7T \J^Ep L p h»=E p [u»+ Xi=2 e^u v \ 

We want to populate the states of these fermions / anti-fermions in this frame according to the 
Fermi-Dirac distribution and calculate the conserved currents in the thermal ensemble. This 
gives 



TilV = E I ^{f q + f- q )Ep[u» + Xd ^u c 



r-OO 

^L{f q + f_ q )Ep[u^ + Xd ^ a u a} 

species 

= EU^u" + p {sT + + CorX + Cn"" 

J "= E J ^ - K + X«^ti„] 

species ^ 
= nil} 1 + Janom 

/•CO JTp 

J 5 = -E y -2^(^+^)K+x d=2 ^o 



(3.2) 



species 

+ Js,anom 



where we have used the relation x d=2 ^ a ' u a X d=2 ^ u \ = 9^" + u^u u . We have collected 
together the deviations from the conventional hydrodynamic expectation under the objects 
with the subscript anom. We get the conventional expressions which could have been naively 
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guessed 



n = 



e = 



s = — 




dE, 



QfqEp 



2vr 



dE p 



2vr 




(3.3) 



where the sum is over every fermionic species with particles and antiparticles counted sepa- 
rately. The anomalous corrections are given by 



At this point, the curious reader might wonder what happens to these relations in the 
case of a Dirac fermion. The massless Dirac fermion is just a left Weyl fermion of charge q 
along with a right Weyl fermion of charge q. In this case, it is evident from the expressions 
above that the conventional terms add up and the anomalous terms cancel out. Hence the 
naive guess turns out to be correct for a massless Dirac fermion. It is not very difficult to 
convince oneself by explicit computation that the naive guess works even for the massive 
Dirac fermion. This then is the first lesson from this exercise : there are transport processes 
in hydrodynamics to which only chiral species contribute . 

Since we will be studying these anomalous contributions in more detail- let us simplify 
them by writing the currents above as 1-forms {q an om, Janom, Js,anom}- Further take a Hodge- 
dual on both sides (which we will denote by an overbar) to remove the e^ u to get 




(3.4) 




(3.5) 
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These are easily calculated from the Gibbs free-energy current 



E 



dE p 

~27T 



■9q X d =2 U 



J a 



dQ a 



djj, 

and q ar 



Js,< 



(3.6) 



+ TJ S ,anom n 



This means we basically have to evaluate only one thermal integral. This can be done by 
pairing up the particle and anti-particle contribution and using the identity 8 



00 dEj 
n "2^ 



p {g q + g- q ) 



-2vr 



2! \2ir) 4! 



21 



This gives 



-2vr 




Xd=2 



species 



(3.7) 



where the sum is performed over the fermionic species with a particle-antiparticle pair con- 
tributing to a single term in the sum. Now we note a crucial fact : the anomalous contribution 
is completely proportional to the £7(1) anomaly coefficient Yl species Xd^Q 2 an d the Lorentz 
anomaly coefficient YLspeties Xd=2 ■ So, we come to the second lesson of this exercise : The 
anomalous transport is completely determined by the anomalies in the system. Hence, we will 
call such a transport as anomaly-induced. 

Let us make this more precise - the anomaly coefficients of a system are neatly summarized 
by a polynomial in gauge field strength F and spacetime curvature 9\. For a collection of 
2d Weyl fermions the anomaly polynomial is given by (see the appendix^ for a review of 
anomaly polynomials for fermions) 



'Panom {F, 9t) — 27T 



F 2 



2!(2vr) 2 



E 



X d =2Q 



species 



4| I / ^ -^d-.! 

'.species 



2d 



where p 1 (9a) is the first Pontryagin class of spacetime curvature defined as 



Pi(K) 



2(2^) : 



ttl 



9h 

' (27T)- 



where we use the notation 



3 See appendix §|X]for a derivation of this/related identities. 



(3.8) 
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Using this, we can write a simple rule to get from the anomaly polynomial to the anomaly- 
induced Gibbs free current 



anom 



(3.9) 



This elegant result tells us that one can just read off the anomaly- induced Gibbs- free current 
of an ideal Weyl gas in 2d directly from its anomaly polynomial. Let us now compare this 



Hence, there is a straightforward algorithm in 2d free fermion theories which takes us from 
the anomaly polynomial of the theory to the polynomial Banom m T an< ^ A 4 which determines 
the anomaly-induced transport. 

One of the main aims of the rest of this article is to generalize this result to arbitrary 
even dimensions. However, in higher dimensions one does not have a powerful principle like 
holomorphy to help us and the anomaly induced transport is hence a more subtle effect to 
derive. So we will spend the next section to formulate a principle which will help us find the 
higher dimensional analogues of the above result. 

4. Adiabaticity of Spectral flow 

One of the crucial lessons one draws from the previous section is that there is a single Gibbs 
free energy current from which energy/charge/entropy currents could be derived. We will 
make the reasonable assumption that this continues to hold true in higher dimensions. In 2d 
we derived the expression for this current in terms of a thermal integral 



where J q should be a 2n — 1 form in d = (2n — 1) + 1 dimensions since it is proportional to 
the Hodge dual of the Gibbs current 1-form Qanom- Taking a Hodge dual on both sides, we 
can remove the overbars and write 





(4.1) 



and a natural generalization of the above expression is 






(4.3) 



Our task is to understand better the physical meaning of Jq and hence formulate some sort 
of an equation for it which then can be solved. 



So what is Jq ? It is clear from the way it occurs in the expression for Gibbs-current 
that Jq is a current whose time-component is just the 1-particle chiral density of states 
participating in anomalous transport. Hence, we will call it the chiral spectral current from 
now on. It tells us how a subset of states related to anomaly flow as the whole fmid(the ideal 
Weyl gas in this case) flows. In 2d, we have 



Hence, two identical states which differ only in their chirality have opposite chiral spectral 
currents. We will assume that this is true in higher dimensions too Jq\d=2n ~ X d = 2n hence 
justifying the adjective chiral. 

In general, we expect the chiral spectral current to be a function of both spacetime co- 
ordinates x and the local rest frame energy E p i.e., Jq = Jq(x,E p ) which emphasizes the 
fact that the chiral spectral current at different energies could be different. Further, the 
subscript q denotes the fact that 1-particle states with different charges can have different 
chiral spectral currents. 

With this understanding, we now turn to the question - how do we determine the chiral 
spectral current Jq for an ideal Weyl gas ? Our crucial tool would be adiabaticity - as was 
argued by authors of [3, 4] (which is a generalization of a d = 4 argument in [1]) the anomaly 
induced transport in arbitrary dimensions is heavily constrained by adiabaticity. At the level 
of states , adiabaticity is basically a statement that the states responsible for anomaly induced 
transport do not get created or annihilated 9 as they move up or down in the energy or in 
case of localized states as they move around in spacetime. This assumption as stated is easily 
formulated in terms of a chiral spectral current - it is just the continuity equation for the 
chiral spectral current in the (x,E p ) space,i.e., 



of states times rate at which their energies are increasing. 

For single particle states with the charge q - there are two forces which lead to a change 
in energy. First is the rest-frame electric force - the work done by the electric force is just 
the electric field times the charge current E^qjq . Second is the pseudo force (since the rest 
frame of the fluid is in general accelerating ) given by —Epa^ where o M is the acceleration of 
the fluid. The work done by it is — a M times the energy current i.e., —a ili E p J q a . Combining 
these together, we get 






(4.5) 



DJ q + —J q E = with J^ = (qE — E p a) A J q 



(4.6) 



It is easily checked that the 2d result J g = ^x d =2 u solves the above equation. 



Except at zero energy due to the anomaly as we will see below 
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5. Chiral spectral current and Anomaly 



A curious reader might wonder how the equation we just derived relates to the adiabaticity 
assumption as it appears in [1, 3, 4] where the discussion was entirely macroscopic with no 
reference to microscopic states. Further, we have not input anywhere the information about 
the anomaly in the above equations. In this section we will show that those macroscopic 
equations could be thought of as arising from the above microscopic equation. Further this 
would also clarify how anomalies are related to the chiral spectral current. 

The flux of density of states is associated with the following energy /charge/entropy cur- 
rents 

poo 

Qanom = / I dEp^JqEpfq 

F JO 

roo 

Janom = ^ / dE p J q qf q 

p JO 

POO 

Js,anom = / / dEp J qH.q 
F JO 

For a general system with anomalies, the statement that these transport processes have 
to be adiabatic is equivalent to the following equation (as shown in [4]) 

Dq a nom + Q, A Qanom E A Janom — T D J g anom + {X \D Janom (^ - ^) 

where 2t[-F] is the anomaly 2n-form in a given theory. This equation above assumes that 
one is working in flat spacetime and hence all gravitational anomalies are turned off. We 
will refer the reader to [4] for a proper derivation of this equation and its consequences for 
a general system. Here we are mainly interested in exploring what it means for the chiral 
spectral current. Substituting the above expressions for the various currents and using the 
following relations 

f = 1 = ®9q_ 

Jl ePify-qp) + I QEp 

Epfq + TU q ~ q/lfq = Q q 
EpDfq + TDUq - QflDfq = 

the adiabaticity equation of [4] assumes an especially simple form 



roo 

E / dE p 

v Jo 



g q DJ q - J 1 



OEp 



+ M 2l = (5.4) 



where as before Jq = (qE — E p a) f\J q . In this form, the relation to the microscopic equation 
that we had derived before is evident. To see this we should integrate by parts - this should 
be done carefully since the boundary contributions do not vanish. We get 



dE p g q 



+^-ewii;=« (5.5) 
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We can now use the continuity equation for the chiral spectral current to set the integrand 
inside the integral to zero. 

Now we turn to the boundary contributions : first the UV contribution - g q falls expo- 
nentially as the rest-frame energy E p — > oo and it is reasonable to assume that the growth 
of J q is slow enough that the contribution from E p = oo is zero 10 . This is consistent with 
the fact that anomalies and their related transport processes are not UV-sensitive despite 
the fact that historically they were discovered in calculations where one should renormalize 
UV-sensitive quantities. 

Now we turn to the IR contribution - there is no good reason for the boundary contri- 
bution from E p = to vanish and in fact we need it to balance the contribution from the 
anomaly. Equating the anomaly to this IR contribution, we get 

a = --2> g (£„ = o) J q E (E p = 0) (5. 6) 

^ F 

This is an interesting equation which is a version of the well-known relation between the 
spectral flow of chiral zero energy states and the anomaly in the system. To see how this 
might work, we will now write the anomaly term also as a sum. For a set of charged Weyl 
fermions in d = (2n — 1) + 1 dimensional flat spacetime, the covariant anomaly is given by 

species 

where the sum runs over every species with one term appearing in the sum for every particle- 
antiparticle couple 11 . 

Let us divide the field-strength F appearing in the above equations into electric and 
magnetic fields in the rest frame of the fluid. The rest frame electric field is E^ = u u F^ v 
which we can think of as a 1-form E. The rest frame magnetic field is a 2-form obtained by 
subtracting the electric part from F, 

B = F — u A E 

where u = u^dx^ is the velocity 1-form. Substituting F = B — E Au into the anomaly above 
we get 

im = E W5 A f (0" ' A ^ (5.8) 

species 

10 In fact, if did grow exponentially with energy the various integrals we have been writing down would 
all be UV divergent and we would have to worry how to make sense out of them. We will assume that this is 
not the case and this assumption can be justified in explicit examples. 

n It docs not matter which among that couple is chosen, since Xd=2nl" +1 ls the same for both - this follows 
from the fact that if the particles charge and chirality is (q,X d =2n) then the antiparticle's charge/chirality is 

(-^(-ir-'x^j • 
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Now use the fact that g q (E p = 0) — g- q (E p = 0) = —q/j, to get 



n—l 



ain_-l£ft(*-O)x«.gA(0 A — jjj (5.9) 

where the sum now runs over all the fermions with particles and antiparticles counted sepa- 
rately. Comparing this with ( |5.6| ), we get for each particle 

J q E (E p = 0) = qE A J q (E p = 0) = g A Xfc2n (g) A 

or 

m = 0) = ^(|f)""A^ (5.10, 

This flow of density of states at zero energy has a direct explanation in terms of Landau level 
physics in the local rest frame of the fluid. We review this connection in more detail in the 
appendix^Q and relate it to the chiral magnetic effect in ideal Weyl gas. Hence, to conclude 
the following equations 

Dj q + —j q E = with 

Jq = {qE — E p a) A J q (5.11) 

implies for an ideal Weyl gas the macroscopic adiabaticity condition in eqn. (|5.2| ). Our next 
task is to solve these equations for an arbitrary fluid flow which we will do in the next section. 



6. Solving the adiabaticity equation 



We now seek a solution of eqn.( 5.11 ) for an arbitrary fluid flow of an ideal Weyl gas. We 
begin by decomposing the exterior derivative of the velocity into its vorticity part and the 
acceleration part 

Du = 2u-uAa (6.1) 

where to is the vorticity 2-form with uj^ u u u = 0. We will now make an ansatz for the solution 
of the form 

_ (2~7r)™ ^ ak ^r~ {n -l-k)\ A U (6 - 2) 

This ansatz is motivated by the solution to the macroscopic adiabaticity equations presented 
in [4] which took the form of a polynomial in 2-forms B and u. The powers of E p are then 
fixed by dimensional analysis. The numerical coefficients are arranged such that the E p = 
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boundary condition fixes ao = 1. We now want to substitute this ansatz into the eqn.(5.2) 
to fix other a k s. Using the following identities [4] 



Du = 2uj + a Au 
D(qB) Au= -qE A 2a; A n 
D(2u>) A« = aA2wAu 



(6.3) 



we get 



DJn 



Xd=2 

(2tt) 



n-l 

E 

fc=0 



qE 

fca fc _i— (k + l)a k a 



A 



(2uE p ) k (qB) 



n—l—k 



k\ 



A 



(n-l - it)! 



A u 



(6.4) 



where we have in addition used the fact that any 2n form made of purely spatial forms B 
and oj is zero. On the other hand 



dE/i 



(2tt)' 



n-l 

E 

fc=0 



qE 

ka k — (fc + l)a fc a 



A (2a;^) fc A J gJ B) 



ri— 1 — k 



in 



l-k)\ 



A u 



(6.5) 



and demanding that the sum of the last two equations should vanish sets a k = a/c-i for all 
k > 1. Along with the boundary condition at E p = which sets ao = 1 this determines 
a*; = 1 for all k. Substituting this into our ansatz, we finally get 



Jo 



Xd=2r 

2n 



qB + 2coE f 
2^ 



n— 1 



A 



n 



1 ! 



(6.6) 



This expression is the central result of the article - it is a formula for how the chiral states of 
a given energy flow when the fluid flows. Since we have not invoked any equations of motion 
for the fluid in our derivation, this is an 'off-shell' solution valid for arbitrary fluid flows - a 
microscopic analogue of the off-shell solution derived in [4]. In the rest of the article we will 
explore various consequences of the above formula. 



7. Anomaly /transport in Ideal Weyl gas of arbitrary dimensions 

We begin by substituting for J q in the expression for Q an0 m 

/■oo 

Ganom — / / AEp^Jq Qq 

p Jo 

- v f°° dEp ( qB+ 2ujEp y 1 u 



(7.1) 



To evaluate this integral, we will again pair the particles and anti-particles together and use 
the fact that if the charge/ chirality of a particle is (q, X d=2 n) then the charge/chirality of the 



-17- 



anti-particle is (-q, (-l) n 1 X d=2 „)- 



Qanom 

= E 



species 



E 

species 



f°° dE p 




Jo 2vr 


5g ( 


f°° dE p 




Jo 2vr 


9q ( 



'qB + 2uE p 

'qB + 2ujE p 
2^r 



n-l 



+ (-i) n "V 



n-l 



+ 9- 



qB - 2uE p 
2^ 



2vr 

n-l' 



n-l 



A 



(n-l)! 



A 



(n-l) 



(7.2) 



To proceed further, it is convenient to employ a formal trick - the trick is to construct a 
generating function which will in one sweep contain in its Taylor expansion Qanom of arbitrary 
even dimensions. To this end let us multiply the above expression by r™ _1 where r is a formal 
parameter and perform a sum over all integers n > 1. We get 

oo 

^ ^ T {Qanom) ( 



71=1 



E E 

n=l species 
species 



1 d=2n 



dE p 
~2T 



9q 



qB + 2w£ p 



dEp 
~2^~ 



9qZ 2 



2tt 



n-l 



+ 5- 



g£ - 2w£ p 
2^ 



n-l 



A 



n — 1 

(n-l)! 



(7.3) 



where we have used Xd t° represent the chirality in the appropriate dimension. 
We evaluate this integral in the Appendix§|A[ The final result is 

CO 

^ ^ T {Qanom) ( 



n=l 



qB 



species 



d=2n 

2ir 



sin £2w/3 



r^ 2 ^ - (i + 



(7.4) 



or 



anomJ d=2n 



n=l 



(2w) 2 ^ 



species 



2ir- 



\2ufi 



sin 52o;/3" 



«- x V 2vr yp 



(7.5) 



We want to now compare this with the form eqn(|ll|) derived in [4] which we reproduce 
for the convenience of the reader 
1 



-[r an(mi [T(2uj),B + n(2uj)] - [$Z nom [T(2cu),B + (,(2u})} 



LU = 



^ nom [T{2u>), B + m(2w)]1 1 

Old Juj=0-> 



(7.6) 



A u 
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By comparing term by term, we get a simple expression 



anom/d=2n 



-2. £ 



species 



IT 

sin 



r n + l 



(7.7) 



where the subscript r n+1 denotes that one needs to Taylor-expand in r and retain the coeffi- 
cient of r n+1 . 

Our aim was to give a prescription to get to danom fr° m the corresponding anomaly 
polynomial. The anomaly polynomial of an Ideal Weyl gas is given by (see appendix [B] for 
explicit expressions in various dimensions) 

(V anom ) d=2n = -27r X d=2n [A(t9K) e^ F 

species 

Now using (see appendix [b] for a derivation) 



sin 



we can write 



(7.8) 



Zanvm = ^anom (F ^ fi, Pl —T 2 , Pk>1 (K) ^ 0) 



(7.9) 



in any dimension for arbitrary collection of free fermions. This is a remarkably elegant 
result which takes the anomaly polynomial and via simple substitutions converts it into the 
polynomial which governs the anomaly-induced transport. 



8. Discussion 



One of our main motivations in undertaking this study was to understand at a more micro- 
scopic level how the anomaly-induced transport comes about. While the transport propor- 
tional to the magnetic field can be understood relatively easily 12 , the microscopic origins of 
the vortical effect were relatively obscure. Having understood now the relevant microscopic 
dynamics as that of spectral flow, one might ask whether one can proceed away from equi- 
librium. In other words, how does this picture of chiral spectral current help us understand 
the non-equilibrium phenomena driven by anomaly. 

While we do not have a complete answer to this question, let us give an idea of what the 
form the answer might take. Ideally to answer non-equilibrium questions, one would like to 
have a Boltzmann-type equation describing an ideal Weyl gas to which interactions can be 
added via collision terms. The Boltzmann equation is intimately tied to the flow of states in 
the semi-classical phase-space : in fact, it is just a conservation-type equation on the phase 
space. Let ( A be the co-ordinates on the phase-space and J^ t be the net flow of states (the 

12 By a Landau level argument which we review in appendix §0 
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total spectral current) in the phase-space. Then the Boltzmann equation takes the schematic 
form 

where the RHS is the collision term which takes into account the scattering from one point 
in the phase-space to the other point. Assuming the total number of states is conserved, we 
can assume dAJtot = so that the Boltzmann equation reduces to 

J t A ot d A f = C[f] 

Hence, we see that a spectral current in the phase space leads us immediately to a Boltzmann 
type equation. To proceed further we need to figure out what the expression for J^ t is ? 
Conventional kinetic theory would suggest that 

x** ~ pt 1 and p^ ~ qp u F tlu 

d d (SI) 

=> JnaivedA ~ j/*^ + ^ F ^ 

which would be identical to the Boltzmann operator one would write down for a Dirac fermion. 
But, as we had argued in this article, in an ideal Weyl gas there is an additional contribution 
to the spectral current over and above the contribution in a Dirac gas. This should lead to the 
modification of the naive Boltzmann operator that we had written above. Unfortunately, the 
spectral current that we have written down in this article has no momentum information - it 
is only a function of {x^,E p } and hence we cannot yet write down the exact modification of 
the Boltzmann operator. It would be nice to work out the momentum-resolved chiral spectral 
current in order to write down this modification . 

We will now give an alternate argument on why we expect such a modification to the 
Boltzmann equation - as is known in the phenomenology of Weyl semi-metals [26, 27] (these 
are 3 + Id phases with Weyl fermionic excitations ), the Weyl fermion in 3 + Id should be 
thought of as a source of Berry flux in the momentum space. The presence of such a Berry- 
flux is known to lead to exactly the kind of modifications of the Boltzmann equations [24, 25] 
that we propose above. It would be nice to clarify the relation between such Berry phase 
related ideas and the ideas presented in this article. 

Let us recap : what we have done is to use adiabaticity as the basic idea which drives the 
anomaly-induced transport in Ideal Weyl gases. This has resulted in a simple rule given by 
eqn.([T.S|) which determines completely the polynomial introduced in [4]. The novel result is 
the way gravitational anomalies seem to enter into the temperature dependence as was noted 
by [11]. While we have reproduced their 4d results, we have done it by an entirely different 
method - they had employed the Ad Kubo formulae derived in [5] whereas our derivation relies 
on adiabaticity. 

Kubo formulae for anomaly-induced transport in higher dimensions involve higher point 
functions and since the Kubo formalism is not well-developed beyond two-point functions, 
it is a priori more difficult to take that route for calculations in higher dimensions. In this 
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article, rather than generalize the Kubo formulae in [5] we have taken an alternate route. But 
having got the answer, let us provide a guess as to how the results of [5] would generalize. 
Consider the small frequency /small momentum limit of the following ra-point function in 
d = 2n dimensions at finite T, fx : 

e Mltt...i»ili2...in-l( T 0ll (k?)T 0t2 (k%) . . . T 0il (kr)Ji l+ M + l)Ji l+ M+2) • • • JiniK) ) 

~ &(T, n) 5 2n (h + k 2 + . . . + k n ) kfk% . . . k^ ] 

then generalizing the results of [5], it is tempting to conjecture that &(T, yu) is the transport 
coefficient associated with (2u) l ~ l A £? n_ * A u term in q an om or equivalently (2w) z A B 11 ^ 1 ^ 1 A it 
term in J an0 m since these two are related by the generalised Onsager type relation [4] 

&Qa.r).nm. $Ja.rt.nm. / Q \ 



5B ff(2w) 

This in particular would mean that for a free theory we expect the transport coefficient 
£l(T,fi) to be related to a Fermi-Dirac integral of the type 13 



It would be nice to show that such a Kubo-type formula as conjectured above does hold 14 , 
since it then opens up the possibility that we can calculate such transport coefficients in 
any field theory. In particular, such a Kubo formula would be of enormous use in checking 
whether the substitution rule in eqn.( |7.9| ) holds beyond free theory. It would be interesting 
to see what happens to this rule as the effective anomaly polynomial is changed via Green- 
Schwarz mechanism. 

There is already some evidence that this rule holds in strongly coupled holographic phases 
[23]. It would be interesting to generalize this by constructing and studying rotating black 
hole solutions in AdS2 n +i with pure/mixed gravitational Chern-Simons terms. This would 
presumably link the proposed modification of Wald entropy [12, 13] in the presence of such 
Chern-Simons terms to the anomaly-induced entropy transport in the dual CFT. 

We do not yet have an intuitive field theory understanding why the higher Pontryagin 
classes do not contribute to the anomaly induced transport that we discuss in this article. 
One way to undertsand this statement better would be to see whether this is true in fluids 
with gravity duals via the calculation we just outlined. One way this might work out is if 
higher Pontryagin classes fall off too fast in the bulk to contribute to boundary transport 
while the first Pontryagin class falls off slowly and asymptotes to —T 2 (2lo) 2 near the AdS 
boundary 15 . 

Another way in which the T 2 coefficients in Ad were calculated was by Vilenkin using 
rotating ensembles [15-20]. This method as Vilenkin used is plagued with technical subtleties 



13 We would like to thank Carlos Hoyos for suggesting such a relation. 

14 This would presumably involve generalising the Kubo formalism to higher point functions[31-36]. 
15 We would like to thank Karl Landsteiner and Mukund Rangamani for discussions regarding this issue. 
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which have to do with trying to define rigid rotation in a relativistic field theory in the 
Minkowski spacetime. But these subtleties can be avoided by studying rotating ensembles on 
a sphere which serves as an infrared regulator. We have done some preliminary calculations 
in Weyl fermions using this method and they seem to agree with the results presented here. 
The details of this calculation will be presented elsewhere. 

Further it would be interesting to generalise all these arguments to spacetimes with non- 
zero curvature. A first step towards such a generalisation would be to first classify the various 
curvature dependent transport processes which are sensitive to anomalies. Since a free theory 
of Weyl fermions is conformal, we can restrict our attention to just the Weyl-covariatised 
curvature tensors [37] and the holographic computations that we had outlined above could be 
easily used to identify the transport processes which are sensitive to anomalies. It would be 
satisfying to understand the chiral spectral current in the strongly interacting systems via 
holography 16 . 
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Appendices 

A. Notes on Fermi-Dirac functions 

Our objective in this section is to evaluate various moments associated with Fermi-Dirac 
16 We would like to thank Jyotirmoy Bhattacharya for discussions regarding this issue 
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distributions. We will start by denning some basic functions which would be useful later on. 
First we want to define the n th poly- logarithm Li n (x) defined via 

-y\k_y 
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the last line being valid for n > 2. The importance of polylogarithm lies in its relation to 
various moments of Fermi-Dirac integrals - it follows from above that 
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/o - /3 fc+r 
We introduce next the n th Bernoulli Polynomials B n (x) - these are polynomials defined 
via the following generating function 

„xt 

(A.3) 
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We now write the identity relating poly-logarithms with Bernoulli polynomials (which itself 
is a special case of the relation between poly-logarithms and Hurwitz zeta function) 
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which implies the following identity for the Fermi-Dirac moments 
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Similarly, we can write 
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Alternately, we can directly write down a generating function for these integrals by using 
the generating function of Bernoulli polynomials. Multiply eqn.(A.5) by r fc , sum over k = 
to oo and then use eqn. (fO) to get 
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which we used in the main text to evaluate the integral in eqn.( [7.4|) . In a similar vein we get 
a generating function for the f q integral 
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These formulae can be used to explicitly calculate Q a 
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The energy /charge/entropy currents can be obtained from these expressions via 

dQ anom 



J anom 

Js ,anom 
Qanom 



(A.16) 
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B. Anomaly Polynomials 



One of the most succinct ways to capture the anomalies in a system is via the anomaly 
polynomial V anom (F, 9t) of gauge field strengths F and the spacetime curvature £R. We will 
not review how these polynomials are calculated or how covariant anomalies are obtained 
from them since these topics are covered well in various textbooks [38-40] and lecture notes 
[41, 42]. We will mainly state the results relevant for Weyl fermions in arbitrary dimensions 
to save the reader the effort of converting the standard results into our notation. 

We will begin by reviewing various forms relevant for dealing with gravitational anomalies. 
Let 9\ a b be the curvature 2-forms of the spacetime with 

1 



Kab = Rabcd dx c A dx d 



(B.l) 



In d = 2n dimensions we can think of 9\ a b as a real 2n x 2n antisymmetric matrix of 
2-forms. At a give point in the manifold, we can diagonalize it with the diagonal entries being 
2-forms 

( +ir 1 \ 



ab 







ir. 



\ 



+ir n 






-ir n J 



(B.2) 



where rj=\ n are real 2-forms. Polynomials in these 2-forms can be used to construct various 
other useful forms. The rest of this section basically consists of various such polynomials, 
relations between them and their generating functions etc. We start with the most basic form 



Ol 



a 2 



. . . A % 



l + (-l) fc 



3 
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As is 



so is a is a 2fc-form which is a /e^-degree polynomial in the curvature 2-forms 
evident from the expression above, it is non-zero only when k is even. 

The next form which is we will introduce is called the fe^-Pontryagin class denoted by 
p k (yt) which is 4/c-form using a specific 2/c-th degree polynomial in the curvature 2-forms. It 
is defined via the relation 
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n 1 =e^(") 
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which can be inverted to get 
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We will express all the other polynomials in the basis of either 9*&s or p k (p$). 

An object that plays an important role in the anomalies generated by Weyl fermions is 
the Dirac genus (or A-roof genus) (9t) which is defined via 



A (r*K) = det 1 / 2 
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which gives 
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The £;* A -Dirac genus A k (9\) is hence a 4fc-form using a specific 2fc-th degree polynomial in 
the curvature 2-forms. 

We now want to calculate the Dirac genus in the special case where the only non-zero 
Pontryagin class is p 1 (iEH) with p k>1 (d\) = 0. It is clear that in this case A k (d\) ~ pf(9t) but 
we want to calculate the numerical coefficient exactly. To do this, we use the following trick 
: we set the 2-forms r 3 >i = with a non-zero r\. Then 
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The roof-genus is easily evaluated in this case as 



1.21- . 
2 2tt J 

'1 T 



sinh lf^ r iJ sinh(§v^(^)) sin (5^/^(91) 



So this gives us the required answer 
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The anomaly polynomial for a bunch of Weyl fermions in a d = (2n — 1) + 1 dimensional 
spacetime is given by 

CPano m ) d=2n = -27T £ Xd . 2n [i (rJK) e^]^ 
species 

where we have assumed that there is a single ?7(1) symmetry under which the charge of the 
Weyl fermion is denoted by the letter q. The sum is over the species which means that each 
particle/anti-particle pair contributes one term to the sum (the answer does not depend on 
whether we take the chirality /charge of a particle or the anti-particle ). We now proceed to 
write down the explicit expressions for V a nom by using the formula above. 
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6!V^ + 4! ^> 
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1 



+- 



2! V2vry 

where explicit expressions for Aj(!SH) are given by eqn.(| 



yl 2 (<K) + A 3 (<K) 



C. Chiral magnetic effect in arbitrary dimensions 

In this appendix, our main aim is to understand in more physical terms how the boundary 
condition for the chiral spectral current arises. In particular, we would like to understand the 
rate at which the zero modes are injected into the fluid in terms of Landau level physics. 
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Let us consider a Weyl fermion in (2n — 1) + 1 dimensions. As before we have a finite 
temperature T, a chemical potential /U and a constant velocity u^. We will again work in the 
rest frame of . The main aim of this section is to understand the transport that arises when 
you turn on a magnetic field in this rest frame. To do this, let us divide the 2n — 1 spatial 
directions into a direction x\ and n — 1 planes where the k th plane is the (xi 2k , Xi 2k+1 ) plane 
where k = 1, . . . ,n — 1. We by consider the system with a uniform magnetic field strength 
turned on in all these n — 1 spatial planes b k = Bi 2k i 2k+1 for k = 1, . . . , n — 1 leaving out 
one spatial direction x\. 

We can solve the Weyl equation in this magnetic background and the solutions are the 
well-known Landau levels with their Landau degeneracies. The energy levels are determined 
by the charge q, the momentum along the x\ direction p\ and the spin along k th plane 
Si 2k i 2k+1 = |<r& (the allowed values are a k = ±1 ) and the Landau level number in the k th 
plane being a non-negative integer n k - 



n-l 



k=l 



\'l ! >i,- ( n k + - J - ^qb k a k 



(C.l) 



where we see the free-particle dispersion along x\ direction, the harmonic oscillator like spec- 
trum of the covariant Laplacian along with the Zeeman splitting due to the magnetic moment. 
The Landau degeneracy in momentum between pi and p\ + dpi (for a fixed Landau level num- 
ber n k and spins a k ) is 



^Landau = ^ ]J ^ (C.2) 
k=l 

We notice that these are effectively a collection of 1 + 1 dimensional free Dirac /Weyl Fermions 
with different masses. So one can repeat the analysis in section §|3] and we conclude that there 
should be a current in this rest frame (along x\ direction) which just depends on the chiral 
states. 

To get a chiral state one has to set the effective 2d mass appearing above to zero. : it 
is easy to see that the 2d- Weyl states are obtained only when n k = 0, a k = sign (qbk)- These 
states have E 2 = p\ = E 2 and a degeneracy 

dpi W \qb k \ _ dEp ( yr qb k 
2tt W 2tt ~ 2vr I H 2vr 

k=l \k=l 

If wg denote the 1 -f~ 1-d.irnensional chircility a-s X.d—2 ^-^^ir totcil contribution to Qanom is 
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where {e\)^ is the unit vector along X\ direction which ensures that the current is along the 
x\ direction. We recognize the appearance of the d = 2n chirality 

n— 1 n— 1 

x d=2n = x d=2 n oh = x d = 2 n sign ( qb ^ 

k=l k=l 

and this equation fixes the 2d chirality x d=2 °f the zero Landau-level in terms of the chirality 
X d=2n °f the original Weyl fermion. This means we get one 2d-Weyl fermion of definite chirality 
for each Weyl fermion in higher dimensions. 

As before, it is useful to take a Hodge dual and write this in terms of forms. Let 
be the field strength which we can think of as a 2-form F. The rest frame electric field is 

= u u F^ v which we can think of as a 1-form E. The rest frame magnetic field is a 2-form 
obtained by subtracting the electric part from F, 

B = F — u A E 

where u = u^dx 11 is the velocity 1-form. In terms of these quantities, we finally get the total 
contribution from all fermions as 

= ^ f 00 dEp (<lE\ n ~ l u 

{Ganom )u=consta,nt — ~2^~^ q ^= 2 ™ J ^ ( n _ fj] (C-4) 

where the sum is performed over all the fermionic particles counting the particle and the 
anti-particle separately. We have denoted in the subscript our assumption that the velocity 
is constant. This gives the rate of injection of chiral zero modes by the magnetic field as 

Certain comments are in order - when we turn on a mild magnetic field over a fluid, 
the gap between the Landau levels (and the Zeeman split) are small and hence one expects 
that any transport will get contribution from many Landau levels not just the lowest Landau 
level. This expectation is in fact correct in the case of conventional quantities like pressure 
or energy density - but remarkably enough the anomaly-induced contribution above does not 
get any contribution from the higher Landau levels even when the magnetic field is small. 
This is related to the statement that only chiral fermions contribute to the sum above . As 
a consistency check, one can easily check that Dirac fermions do not contribute to this sum 
( this follows since in a Dirac species for a given value of q , X d=2n takes both the values ±1 
thus canceling out in the above sum ). 



D. Notation 

We work in the ( !-+•••) signature. The dimensions of the spacetime in which the fluid 

lives is denoted by d = 2n. The Greek indices /x, v = 0, 1, . . . , d — 1 are used as space-time 
indices. 
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We denote Hodge-duals by an overbar - for example, J is the 2n-l form Hodge-dual to 
the 1-form J„. We mostly just use the Hodge-duality between 1-forms and 2n-l forms and 
our conventions are completely defined by the following statement- given any 2n — 1 form V 
hodge-dual to and a 1-form A^, we have 

DV = (D,Vn Vohn 
A AV = —V A A = A^V No\ 2n 

Given a 0-form a its Hodge-dual 2n-form is simply a = a Vol2 n . 

We have included a table with other useful parameters used in the text. In the table 
the relevant equations are denoted by their respective equation numbers appearing inside 
parentheses. 
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Table of Notation 


Symbol 


Definition 


Symbol 


Definition 


e 


Energy density 


P 


Pressure 


n 


Charge density 


s 


Entropy density 




Chemical potential 


T 


Temperature 


a 


1/T 




2n dimensional chirality 


q 


Fermion charge 


E p 


Fermion energy 


9q 


-Am 

PL J 


fq 


\ e l3{E v -qa) + 1 1 - 1 
L J 




Curvature 2-forms(B.l) 




See 

\r V 


Tanom 


Anomaly polynomial 




Pontryagin class flB.4|) 


'Mt) 

u anom 


See (1.3) 


A, (fK) 

k \ J 


A-roof genus QB.8|) 




Fluid velocity, 
1-form 


CLi,, CL 


Acceleration field 
(u.D)u a , 1-form 




Spacetime metric 


p 


G,,,, -\- ?/,,?/>/ 




Shear strain rate 




Fluid vorticity, 2-form 




Energy-momentum 
tensor of the fluid 




Ohrir0'p nirrpnts 
with anomalies 


°s 


Entropy current 


D 


Exterior derivative 


yanom 


Anomaly- induced 
Gibbs current 


Qanom 


Hodge-dual of Qanom 

In — 1 form 


a 

Qanom 


Anomaly- induced 
heat current 


Qanom 


Hodge-dual of qanom 
2n — 1 form 


TV- 
's anom 


Anomaly- induced 
Charge current 


J anom 


Hodge-dual of Janom 

2n — 1 form 


b,anom 


Anomaly- induced 
Entropy current 




Hodge-dual of J a 

ii^^ 6 ^ uuui " S,anom 

2n — 1 form 


F F 


non-dynamical gauge 
field strength, 2-form 


E^,E 


Rest frame electric 
field F au u v ', 1-form 


Bfj. u , B 


Rest frame magnetic 
fields F — u A E 


21 


Anomaly 
DJ = 21 


Jq 


Chiral Spectral current 


J* 


Hodge-dual of Jq 


jE 


See eqn([f .5[) 


n-E 


Hodge-dual of X E 
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